Phase diagram of iron-pnictides if doping acts as a disorder 
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We obtain and analyze the phase diagram of doped iron-pnictides under the assumption that 
doping adds non-magnetic impurities to the system but does not change the densities of carriers. 
We show that the phase diagram is quite similar to the one obtained under the opposite rigid band 
assumption. In both cases, there is a phase where superconductivity and antiferromagnetism 
co-exist. We evaluate the jump of the specific heat, AC, at the superconducting Tc across the phase 
diagram and show that AC/Tc is non-monotonic, with the maximum at the onset of the co-existence 
phase. Our results are in quantitative agreement with experiments on some iron-pnictides. 



I. INTRODUCTION 

How chemical doping of iron-pnictides affects their 
electronic structure is not fully understood yet and is 
subject of debates. In most studies it is assumed that 
doping does not affect the rigid band picture and only 
changes the densities of holes and electronsi. An al- 
ternative scenario^ is that doping does not affect the 
carrier density but rather introduces non-magnetic im- 
purities and hence increases disorder. Angle-resolved 
photoemission (ARPES) experiments on 122 materials 
Ba(Fei_a;Cox)2As2 and Bai_a;Ka;Fe2As2 are usually in- 
terpreted in favor of the rigid band scenario. Within 
this scenario, if magnetism prevails at zero doping, the 
system moves from a spin-density-wave (SDW) phase to 
a superconducting (SC) state, and for some model pa- 
rameters there is a mixed phase, where SDW and SC 
orders co-exist.—"— Recent ARPES experiments on Ru- 
doped BaFe2As2, however, found^ii that substitution of 
Fe with Ru practically does not change the Fermi sur- 
face (FS), yet the phase diagram is quite similar to that 
in other doped 122 materials: as Ru concentration in- 
creases, the system moves from an SDW phase to an SC 
phase. In between, there is a region where where both 
SDW and SC orders co-exist, although microscopic co- 
existence (as opposed to phase separation) has not been 
experimentally proven yetj^ Because FS geometry does 
not change, it seems natural to assume that the changes 
in the phase diagram caused by Ru-substitution are pre- 
dominantly due to dilution and disorder associated with 
it. We also note that disorder may be introduced directly 
to pnictide materials by irradiationj^iiS 

In the present paper, we address the issue of what is 
the phase diagram of a doped 122 Fe-pnictide if doping 
does not affect the carrier density but rather introduces 
non-magnetic impurities. We show that the phase dia- 
gram is actually the same as in the rigid band scenario. 
Namely, as doping increases, first an SDW phase becomes 
a mixed phase, then the system becomes a pure SC, 
and at even larger dopings SC is destroyed by disorder. 
This result may look somewhat counter-intuitive because 
non-magnetic impurities are pair-breaking for an SC. 
It turns out, however, that impurities damage SDW order 
stronger than they damage SC because both intra and 



inter-band impurity scattering is destructive for SDW^ 
while only inter-band scattering is pair-breaking for an 
SCii^rii Because of this disparity, the actual mag- 
netic Ts becomes smaller than the superconducting Tc 
when the density of impurities exceeds a certain thresh- 
old, even for the undopcd case > Tc- There is no 
a'priori guarantee that a mixed state emerges near the 
point where Tg = Tc, i.e. a first order transition from an 
SDW to a SC is another option. Our calculation shows 
that the mixed state does appear, see Fig. [T^. For such 
a phase diagram to emerge, the magnetic SDW Ts.q for 
undoped material should not be too strong compared to 
Tcfi, see below. If Tsfi/Tc.o is too large, Ts remains higher 
than Tc down to Tc —> 0, even though T, decreases faster. 

There is another reason to analyze the phase diagram 
assuming that doping introduces disorder. The measure- 
ments of the specific heat jump AC at Tc across the 
phase diagram have dcmonstratedi^"— that AC/Tc is 
non-monotonic and has a maximum at optimal doping 
that almost coincides with the onset of the co-existence 
phase. The slopes of AC/Tc are similar, although not 
exactly identical, upon deviations from optimal doping 
into both directions. This similarity raised speculations 
that the behavior of AC/Tc in underdoped and overdoped 
regimes may be related. Within the rigid band model, 
AC/Tc has a peak at the onset of a mixed phase and 
rapidly decreases at lower dopingii^ However, the reduc- 
tion of AC/Tc at higher doping cannot be straightfor- 
wardly explained within the rigid band model. 

In the disorder model, the behavior of AC/Tc across 
the whole phase diagram is determined by a single pa- 
rameter, the density of impurities nimp, and the forms 
of AC/Tc in the underdoped and overdoped regimes are 
related. We find that in the disordered model, AC/Tc in- 
deed decreases on both sides of optimal doping, as shown 
in Figs. [1)3 and [21 The specific heat is discontinuous at 
the onset of the mixed phase within the mean field ap- 
proximation, but becomes rounded once fluctuations are 
taken into account. The decrease of AC/Tc away from 
the maximum shows rather similar, although not iden- 
tical, behavior in under- and over-doped regimes, with 
roughly quadratic dependence on the transition temper- 
ature Tc, see Fig. ^p. This behavior is in quantitative 
agreement with experimentsi^"— 

The fact that the phase diagram and the behavior of 
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AC/Tc are similar in the rigid band and the disorder 
models is encouraging, since the two models are comple- 
mentary to each other. In general, a chemical doping acts 
in both ways: (1) doping introduces some extra carriers 
and (2) increases impurity density. The relative magni- 
tude of the two effects depends on materials. We argue in 
this regard that quite similar behavior observed in Ru, 
Co, and K - doped BaFc2As?i^"— is not a coincidence 
but rather a quite generic feature of iron-pnictides. 

The paper is organized as follows. In the next section 
we discuss the model and introduce SDW and SC order 
parameters and describe the formalism used for calcula- 
tions. In Sec. mil we analyze the phase diagram as a 
function of impurity concentration, by solving linearized 
gap equation for one order parameter, SDW or SC, when 
the second parameter is either absent or present. Sec- 
tion |IV] presents calculations of the superconducting or- 
der parameter near the transition to a superconducting 
state. In Sec. |Vl we consider specific heat jump at the 
onset of superconductivity. We provide our conclusions 
in Sec. Ell 



II. THE MODEL 



erators and can be written as 



(1) 



where ^p,Q — (^p^a, ^-p^a^ /p,Q> f^p,a) and *p,Q is a 
conjugated column. The Hamiltonian matrix Hp^a.p can 
be written in the formi 

H = Hq + Hmf, Hq = -^Tsps&o 

Hnif = -Af2P3(T2 + T3Pi{Ma-). 



(2) 



Here, the Pauli matrices t^, pi and (t, are defined in 
the Gorkov-Nambu, band, and spin spaces, respectively, 
where i = 0, 1, 2, 3 and matrices with i = are unit ma- 
trices. A fermion Green's function G{oJn,p) is defined as 
a solution to 

(^iujn - Hp -'Ej G{LLin,p) = 1, (3a) 
and the conjugated equation is 

G{^n,p)(i^n-Hp-t) ^i. (3b) 



A. General formulation 

Our goal is to demonstrate that the phase diagram 
remains the same if we associate doping with disorder 
rather than with the changes to the FS in the rigid band 
picture. We adopt the same minimal model that was 
used in earlier works within the rigid band approach.-i^ 
Namely, we consider a two band metal with cylindrical 
FSs for electron and hole-type excitations. The cylindri- 
cal FSs have circular cross-sections of equal radii centered 
at (0,0) with a hole-like dispersion and Q = (0,7r) with 
an electron-like dispersion. The free fermion part of the 
Hamiltonian in this case of perfect nesting is represented 

by 



Here S is the self energy for scattering off disorder and 
ujn = 2'KTm{n + 1/2) with integer n are Matsubara fre- 
quencies. 

We describe disorder scattering within the Born ap- 
proximation and assume that the Born scattering ampli- 
tude U (q) is characterized by a constant Uq for scattering 
within the same band and J7,r for scattering between the 
two bands.— ^— In this approximation, the self-energy is 



dp 
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where we introduced disorder scattering rates 
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where operators c annihilate hole-like fermions near (0, 0) 
and operators / annihilate electron-like fermions near Q. 
The fermionic dispersion is given by ^(p) = p^/2m — /i, 
and the momentum p of electron excitations is measured 
as a deviation from Q, p ^ p — Q. 

We consider an effective low-energy theory with the 
high-energy cutoff A and angle-independent interactions 
in the SDW channel and in the SC channel with the 
couplings Asdw and Asc'^'^fli"— We treat these interactions 
within a mean field approximation, by introducing SC 
and SDW order parameters, A and M, respectively, and 
decomposing the four-fermion interactions into effective 
quadratic terms with A and M in the prefactors. The 
full mean-field Hamiltonian is quadratic in fermionic op- 



Fo characterizes the rate of electron collisions with impu- 
rities in which the electron remains in its original band, 
while F^ is the rate of collisions resulting in electron 
transfer between the two bands. Np in Eq. ([5]) is the 
total quasiparticle density of states (DoS) at the Fermi 
energy ( the DoS per spin per band is Np/A). We assume 
that only the impurity density riimp changes with doping, 
i.e., the ratio F^/Fo is doping independent. 

The two mean-field parameters A and M are obtained 
self-consistently via the matrix Green's function as 



A_ 



T 



T.jj^^^{^i-r.,p)fUpO + P3)a^} 



(6) 
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and 
M 



sdw 



lE//(^Tr{GK,p)(fo + r3)p+<.}, 

where = {Ai +iAi)l2 for A p,f,a. 

For the pure SDW and the pure SC state 

in the absence of disorder, the solution of the 
hnearized gap equations yield transition tempera- 
tures Ts,o = 1.13Acxp(-2/(A^FAsdw)) and Tc,o ^ 
1.13Aexp(-2/(A^FAsc))- We consider T^^ > Tc,o, so that 
without disorder the SDW phase develops at a higher 
temperature. 



B. Eilenberger equation 

To treat superconductivity and magnetism in the pres- 
ence of disorder, it is convenient to introduce the Eilen- 
berger's Green function 



4i 



dp 



(8) 



which appears both in the self-consistency equations, 
Eqs. (|6]) and (O, and in the expression for the impu- 
rity self-energy, Eq. 1^. In particular, the impurity self 
energy is 

S = -iToTops&ogTspQaQ - ir^-ro(-ip2)o-o^T3/5i<To- (9) 



To derive the equation for Q, we multiply Eq. (jSaj) by 
73 Pa from left and subtract Eq. (|3b|, multiplied by r^p^ 
from right. We then multiply the resulting equation by 
fspa from left again. As a result, the Ho{p) term falls out. 
We integrate the resulting equation over p and obtain the 
equation for ^(w„) in the form of a commutator: 



i'^nTaPs^Q - iH„,i + i;)f3/53(To; = 0. (10) 



This equation is the Eilenberger equationj^i^ obtained 
for a two-band metal with homogeneous in space SDW 
and SC order parameters. The Eilenberger equation 
is consistent with the normalization relations for Q: 
Tr^(a;ii) = and (7(w„)C?(a;„) = 1. 

Without loss of generality, we direct M along z-axis 
and parametrize the matrix t?(w„) by the three functions 
guj„ , fuj„ and S'„„ as 

G{ujn) = guj^hpsaa + /w„Ti/Oo(-i(5-2) + S'aj„ to(-«P2)ct3. 

(11) 

The function g^^^ is the normal component of the Eilen- 
berger Green's function, while the functions S^^ and f^^ 
are the two anomalous components, associated with the 
SDW and SC orders, respectively. 

For the above parametrization of C7(cli„). Eq. (ITT|) . the 
normalization condition QQ ^ I reduces to 



9u. 



f 



(12) 



and the commutation relation Eq. (jlO[) gives 

i/^gu,„ = fu„i^n + 2r^g^„), (13a) 
iMg^^ = (w„ + 2rtg^J. (13b) 

where Ft = Fq 4- Fjr. 

The self-consistency equations for SDW and SC order 
parameters, Eqs. ^ and ([7]), can be rewritten in terms of 
anomalous SDW and SC components of the Eilcnberger's 
Green function as 



2M 
2A 



= -i 2ttT ^ Suj„, 
A 

= -i 2ttT V L . 



(14a) 
(14b) 



III. PHASE DIAGRAM 

We first consider pure SDW and SC states. For a pure 
SDW state we set A = and = in Eqs. ^ and 

linearize Eq. (|13bp in M and find from Eq. 
that the SDW transition temperature evolves with 
doping as 



2-kT. 



sdw 



n>0 ^ ' 



2Tt 



(15) 



This equation can be rewritten in terms of the transition 
temperature Tsfl to SDW phase at Ft = as 



T, ^ I 2 ^T, 



(16) 



where ij^ix) is the digamma function.— 

For a pure SC state we set M = and S't^„ = in Eqs. 
(fT2)) and (jl3bp and linearize Eq. ([T3a| in A. We obtain 
from Eq. (jl4bp 



NfXs 



= 2ttT, 



A/27rTe 

^ 7rTc(2n+l) 



2F„ 



(17) 



Re-expressing the result in terms of the superconducting 
transition temperature T^^q in a clean system and without 
SDW, we re-write Eq. dH]) as 



In 



-nTr 



(18) 



which is similar to the equation for Tc in conventional 
s— wave superconductors with magnetic impurities^l and 
in unconventional d— wave superconductors with poten- 
tial impurities;^"— Note that only inter-band scattering 
Fjr, is pair-breaking for SC. 

Even if Tg^o > ^c.Oj decreases faster than Tc with 
increasing riimp, and at certain doping the two transition 
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FIG. 1: (Color online) Upper panel: The phase diagram as 
a function of disorder, measured in units ro/27rTc,o 
for on-site disorder (F^ = Fo), and Tsfi/Tc,o = 1.7. The four 
transition lines terminate at a tetra-critical point P, where 
normal, pure SDW, pure SC, and mixed phase meet. The 
shaded region represents the mixed phase. Lower panel: Spe- 
cific heat jump AC/Tc as a function of doping. Solid curves 
represent the mean field result and the dashed curve illus- 
trates the effect of thermodynamic fluctuations beyond the 
mean field theory. 



temperatures may cross. We denote this temperature 
as Tp. The condition that Tp exists, i.e. that Tc and Ts 
cross before Tc — >■ 0, sets the limits on the ratios Tgfi/Tcft 
and T^/Tq. We find that and cross if Tc,o/T„o > 
1/(1 + Fo/FTr). For on-cite disorder potential F^ = Fq, 
and Tp > exists, i.e. SC phase exists, if Tc^o/T's.o > 1/2. 
For longer-range impurity potentials, Fq > F^, and the 
SC phase develops even for smaller Tcfi/Ts_o, see Figs. 
and[3^. 

To obtain the superconducting transition temperature 
Tc(M) in the presence of pre-existing magnetism one has 
to solve the linearized equation for A at a finite M . Now 
g^j^ depends on M (i.e., Qu)^ = gi^^{M)), and we have 
from Eq. (|14bp 



where 



a(r,(M)) 



w„>0 



NpX, 



9u,„ (M) 



2F,5„„(Af)- 



(19) 



(20) 



The temperature dependence in the r.h.s. of Eq. ([^0)1 
is via u)n = TTT{2n + 1) and also via g^^^^M) because 
M depends on temperature. The summation in a(T) 



has a logarithmic dependence on the high-energy cut- 
off A. This dependence can be eliminated in favor of 
the transition temperature Tc at M ~ 0. Subtracting 
Eq. p7)) from Eq. ([T^ . we obtain after a simple algebra 
an equation on Tc{M) in the form 



C{Tc{M),Tc,T^) = 



2^Tc{M)L0r.[guAM) ~ I] 

(w„ + 2F,)(cj„ + 2F,5^„(A/))' 

(21) 



where 

(22) 

We calculate g^^{M) as a function of temperature at 
a given impurity concentration. For this purpose, we 
express Si^^{M) in terms of gi^^{M) using Eq. (|13bp . 



iMg^SM) 
Wn + 2Ft5^„(M)' 



(23) 



substitute the result into Eq. (fT2|) with fuj„ = 0, and 
obtain the fourth-order algebraic equation for g^^ (M) 
as a function of M: 



M^gl(M) 



{L0n + 2Ttg^^{M)y 



= 1. 



(24) 



We solve this equation, obtain g^^{M), substitute the 
result back into Eqs. (|23|) and ()14ap . utilize the definition 
of Ts, and obtain the non-linear equation for M ~ M{T) 
in the form 



£(r,T„Ft) = 2^r 



LUn[gu,SM) ~1] 



^>0 



(L^„ + 2F0(w„ + 2Ft(7<,„(M))' 

(25) 

where gui^{M) is a solution of Eq. ([24]), one has to choose 
the branch with g^^{AI = 0) = 1. Solving we obtain 
M(r), and hence gui^{T). Substituting the result into 
([T9| and ([20]) we obtain the superconducting transition 
temperature Tc{M) in the mixed phase as a function of 
doping. 

We numerically evaluate Tc{M) at different dopings 
in the mixed phase and plot the result in Figs. [T^ and 
|3^. As the doping decreases from its optimal value, M 
increases at a given temperature T, and Tc{M) rapidly 
drops. This is expected since SDW and SC order param- 
eters compete with each other. At M — >■ 0, g^^{M) — 1 
and Eq. yields Tc{M) T^ as expected. 

A similar calculation of the SDW transition tempera- 
ture from the preexisting SC phase, Ts{A), shows that 
Ts{A) decreases as A increases, due to the same kind 
of competition. Furthermore, the Ts{A) curve actually 
bends toward smaller dopings, see Figs. [1^ and [3^, so 
that with decreasing temperature the system moves from 
a pure SDW magnet to a pure superconductor through 
a mixed phase. The bending of the T's(A) curve is in 
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agreement with the general analysis in Ref. HT]. The four 
curves T^, Tc{M), Tg, and Ts{A) meet at the tetracrit- 
ical point P, as shown in Figs. [T^ and The corre- 
sponding temperature Tp is the highest superconducting 
transition temperature. We also see from numerics that, 
despite bending, the curve Ts(A) is always located to the 
right of the curve Tc{M), i.e. if one increases disorder at 
a given T or decreases T at a given disorder, the sys- 
tem with the SDW order first becomes unstable toward 
an intermediate mixed phase where SDW and SC orders 
co-exist, and only then SDW order disappears. 

The intermediate mixed phase was earlier found in the 
rigid band modeli^i^ However, in that model it only ap- 
pears at a finite ellipticity of electron pockets, while for 
circular hole and electron FSs doping gives rise to a first 
order transition between pure SDW and pure SC phases. 
In the disorder model, the mixed phase appears already 
for circular hole and electron pockets and by continuity 
should also exists when electron pockets have weak el- 
lipticity. We, however, did not analyze the whole range 
of ellipticities and therefore cannot exclude a possibility 
of a first order transition for strongly elliptical electron 
FSs. 




FIG. 2: (Color online) The specific heat jump AC/Tc for r,r = 
To, and Ts,o = l-TTcO, as a function of Tc/Tp, where Tp is 
the temperature of the tetra-critical point. Left panel: linear 
scale, right panel: log-log scale. The solid line represents the 
specific heat jump AC/Tc at Tc{M) in the mixed phase in the 
underdoped regime, and the dashed line represents AC/Tc in 
the overdoped regime, where SDW order is absent. A thin 
dash fine represents a quadratic dependence AC/Tc oc T^. 



IV. SUPERCONDUCTING ORDER 
PARAMETER NEAR Tc 

We verified that the mixed phase does indeed exist in 
the disorder model with circular FSs by expanding in 
Eq. (|13ap to order and solving the equation for A 
in the presence of M at a temperature slightly below 
Tc{M). The expansion yields, quite generally: 



a{T) - /3A2 



Xsc Nf 



(26) 



where a(T) is introduced in Eq. 1^ and P = I3{T^{M)) 
is given by Eq. ^\ below. Near Tc{M), we have a{T) = 
a{Tc{M)) + a'{T^{M)) (T - T^{M)) and a{Tc{M)), see 
Eq. p9p . On general grounds, a'{Tc{M)) must be neg- 
ative for a SC phase to develop as T decreases, and we 
indeed show below that a'(Tc(M)) < 0. The type of the 
transition is, however, determined by the sign of /3. The 
mixed phase exists if /3 > because then A gradually 
grows as T decreases. If /3 < 0, A changes discontin- 
uously around Tc(M) and the SDW and SC phases are 
separated by the first-order transitioni^ii 

The coefficient a{T) can be rewritten in the form of 
Eqs. (ini) and 



a{T) 



XscNf 



Y{T), a{T) 



dY(T) 
dT 



where 



Y{T)=C{T,T,,T^) 

, ^ 27rTw„(l-5„„(Af)) 



j-^^ [ojn + 2r^][w„ + 2r^5„„(A/)] ■ 



(27) 



(28) 



For i\/ = and in the clean limit a'{Tc) = -l/Tc- We 
verified numerically that a' {Tc{M)) remains negative at 
M 7^ and in the presence of disorder, as expected. 

Calculations of /3 require more care as one has to com- 
bine terms coming from the apearance of non-zero fu)^ 
in Eq. and from the expansion of the SDW order 
parameter M(A) to order A^ as M(A) = A/ + 5M'^'^\ 
where JM^^^ oc A^. Similarly, we introduce g^j^ = 
g„„ {M) + dg^J} and S^^ = S'„„ (M) + SSj^f]. Substituting 
g^^ and into Eqs. ((T2)) and l|13bp . we obtain equa- 
tions for Sq^i^} and SsHf]: 



[un + 2Ttg^„ (A/))2 + ujn + 2Ttg^„ (M) 



(29) 



where 



(1) 



tAff^jM) 
+ 2r,5^„(A/) 



(30) 



and guj,SM) is defined by Eq. Solving Eq. ^ we 

obtain 

^^^(2) ^ _i 9^AM){i^n + 2rt5^„(M))2 



X A 



2(w„ + 2rt5„„(M))3 + M2w„ 

u}n + 2rtguj„iM) 2Af(5A/(2) 



(w„ + 2r,g„„ (Af)2 w„ + 2Ttguj„ (A/) 



(31) 
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and 

c<.(2) ^ ^ g^„(Af)K + 2rtg^„(M))^ 
2{iOn + 2Ttg^jM)r + M^ujr^ 



K + 2r^g^„(M))2(w„ + 2rtg^„(A/)) 



(32) 



We first evaluate SM^"^^ by substituting Sui„ = 5"^°^ + 
Sj^} into Eq. (|14ap and eliminating the SDW coupling 
constant in favor of the SDW transition temperature , 
see Eq. (fTSj) . We obtain 



w„>0 



+ 2rt 



Substituting C(T,Ts,Tt) 



Eq. ([25]) and Ss!;^,} from Eq. ([32]) we obtain 



(33) 

C{Tc{M),Ts,rt) from 



2A/ S ' 



(34) 



where the coefficients B and C, together with the term 
A which we utilize below, are given by 



A = nTY: 



u:n>0 



K + 2r,g^„ (M)]4Z? 



ij„>0 



(35a) 



(35b) 



^ ^ .9^„(^4")^n(w„ + 2rt5<^„(A/)) 



w„>0 



K + 2r^5<.„ (M)]21? 



= K + 2rtg<,„ (Af )]^ + w„Af 



#2 



(35d) 



Substituting JA/^^) into we obtain (5.g(2) A^. 

We next write f^^, defined by Eq. (|13ap to the third 
order in A 



lAg^^iM) 



iAuJng 



(2) 



uJn + 2r,g„„ (A/) (w„ + 2r^.g„,. (1/))^ 



r, (36) 



substitute this expression into Eq. (|14bp . and obtain 
Eq. (HH) with 



B 



(37) 



where B, and C are given by Eq. ([35]). 

Evaluating these coefficients, we find AB > C^, i.e., (3 
is positive. This confirms our numerical result that the 
phase diagram of the disorder model contains the mixed 
phase where SDW and SC orders co-exist. 

Equation ([25)1 can also be applied to the transition 
from a paramagnetic metal into a pure SC state. In this 
case, M = 0, and hence g^^ = gt^^{M) = 1 and C"^ jB 



term in Eq. ([37]) is absent, i.e. /? = A. Substituting 
g^^^ (M) = 1 into Eqs. pSap and p5d[) we obtain 



1 / °° 



E^/ttTc 



/2 + rj7rTe) 



(38a) 



where Tc is given by Eq. Note that, again, q;'(Tc) < 
and ;3 > 0. 



V. SPECIFIC HEAT JUMP AT THE ONSET OF 
SUPERCONDUCTIVITY 

The specific heat jump at Tc and Tc{M) can be ob- 
tained by evaluating the change in the thermodynamic 
potential dfl imposed by superconductivity^ 



An 



dAi 



■AldAi = - 



NfPA^ 



(39) 



where dX^J/dA and A^ are defined by Eq. For A^ 

we have 



A^ = 1 ( a{T) - 



Xsc Nf 



'■^{Tc{M)-T). (40) 



The change of the specific heat due to the superconduct- 
ing ordering is AC = -Td'^Sn/dT'^. At T = Tc{M), the 
specific heat exhibits the jump given by 



AC ^ 37 [a'{Tc{M))Y 

Tr 



27r2 



/3 



(41) 



where 7 = n'^Np/i is the Sommerfeld coefficient in the 
metallic phase. 

The behavior of AC/Tc as a function of doping- 
induced disorder is shown in Figs. [T}d, [21 and ^jp. To 
evaluate AC/Tc at the transition from a normal metal to 
a superconductor above optimal doping we use Eqs. ([55)1 
for a' and /3 in Eq. (|1T|) . At large doping, when Tc is sig- 
nificantly suppressed and the system enters the regime 
of impurity-induced gapless superconductivity, specific 
heat jump AC/Tc decreases with Tc as AC/Tc oc T^, 
Ref. [33. Away from the gapless regime, the dependence 
of AC /Tc on Tc is more complex and differs from T^^ as 
the dashed lines in Figs. [T]3,[21 and[3]3. For completeness 
we present in Fig. [J] AC/Tc as a function of Tc/Tc.t^ for 
an superconductor, when there is no competing SDW 
instability (i.e., when Asdw = 0) and the Tc line extends 
to Tcfi in the clean limit. 

For the transition from the preexisting SDW state 
into the mixed state below optimal doping, we compute 
AC/Tc, Eq. ([31]), using a' and /3 from Eqs. ([17]) and 
([57]) . In this regime, SDW order strengthen as doping 



FIG. 3: (Color online) a) The phase diagram as a function of 
disorder, ro/27rT'c,o for = 0.2ro and Ts.o/Tc,o = 3. Note 
that the region of the mixed phases gets narrower than the 
region for on-cite impurities, cf. [1] b) The specific heat jump 
as a function of Tc/Tp. The solid line represents the spe- 
cific heat jump AC/Tc for normal-to-SC transition, and the 
dashed line shows AC/Tc in the overdoped region for transi- 
tion between SDW and SC coexistence phases. 



decreases, and SDW correlations suppress both super- 
conducting Tc{M) and AC/Tc- In particular, the rapid 
decrease of AC /Tc at smaller dopings is an indicator that 
fewer quasiparticle states participate in superconducting 
pairing, as the low-energy states are pushed away from 
the FS by SDW order. We see therefore that AC/Tc 
drops at deviations from optimal doping in both over- 
dopped and underdopped regimes. 

It is essential that in the disorder model, the behavior 
of AC /Tc in the underdopcd and overdoped regimes is 
governed by the single parameter Fp oc riimp, assuming 
that the ratio Ftt/Fq is kept constant. The same parame- 
ter Fo also defines T^/Tp, where Tp = T^M 0) is the 
transition temperature at the tetracritical point. One can 
therefore make a direct comparison with experiments by 
plotting AC/Tc above and below optimal doping as the 
function of the same Tc/Tp using AC defined by Eq. (|4T|) 
with a' and /? given by Eqs. ([55)1 above optimal doping, 
and by Eqs. ([TT)) and ([57]) for a finite M below optimal 
doping. Experiments showi^"— that AC /Tc drops faster 
in underdopcd regime, but in log-log plot the data from 
underdopcd and overdoped regimes can be reasonably 
well fitted by a quadratic law AC/Tc oc T"^- 

We plot our AC/Tc in Fig. [5] as functions of Tc/Tp in 
both linear and log- log plots. We see that AC/Tc drops 
faster with decreasing Tc in the underdopcd regime, i.e. 
in the mixed phase. This behavior is consistent with 
experimental data. The T dependence of AC /Tc is not 
exactly T^ , but looks reasonably close to in log-log 
plot (right panel in Fig. [5]), although such plot puts more 
weight on the data at low Tc where the T— dependence 
of AC/Tc is the strongest. 

Finally, we note that AC /Tc is discontinuous at the 
tetra-critical point Tc = T^ = Tp. The discontinuity is 
the manifestation of the discontinuous change in both a' 
and P across Tp. The discontinuity in a' is due to the 



FIG. 4: (Color online) The jump AC /T^ of the specific heat 
at a superconducting Tc without a competing SDW instabil- 
ity (Aadw = 0). AC/Tc is plotted as a function of Tc/Tc.o, 
where Tc.o is the superconducting transition temperature in 
the clean limit. 

fact that the second term in (P?]) is zero for M = 0, but is 
finite when M ^ and contains temperature derivative 
of - 1) oc M2 cx (T - Tp). The discontinuity in /3 is 
due to the feedback term C^ / B in (3, which remains finite 
upon approaching Tp from smaller dopings, but is absent 
in the overdoped regime, where Tg < Tc. The interplay 
between discontinuities in a' and /? in the disorder model 
is such that AC /Tc jumps up at riimp once the system 
enters the mixed phase. 

The discontinuity in AC /Tc at Tp has also been found 
in the rigid band modeljiS In that model, however, the 
magnitude and the sign of the jump in AC /Tc depend 
on the FS geometry and AC /Tc may actually drop down 
upon entering into the mixed phase. We also emphasize 
that discontinuity in AC /Tc only holds within the mean- 
field theory and gets rounded up and transforms into a 
maximum once we include fluctuations because then the 
thermodynamic average (Af^) is non-zero on both sides 
of the tetra-critical point. This behavior of AC/Tc in 
the presence of fluctuations is schematically illustrated 
in Fig. [1)3 by the dashed line. 



VI. CONCLUSIONS 

In this paper we obtained the phase diagram of doped 
Fe-pnictides and the specific heat jump AC at the onset 
of superconductivity across the phase diagram under the 
assumption that doping introduces disorder but does not 
affect the band structure. The phase diagram is quite 
similar to the one obtained in the rigid band scenario 
and contains SDW and SC phases and the region where 
SDW and SC orders co-exist. The ratio AC /Tc, which is 
a constant in a BCS superconductor, is non-monotonic 
across the phase diagram - it has a maximum at the 
tetra-critical point at the onset of the mixed phase and 
drops at both larger and smaller dopings. The behavior 
at large and small dopings is described in terms of the 
single parameter: the impurity density njmp. The non- 
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monotonic behavior of /^C /Tc in the underdoped regime 
also holds in the rigid band model^^ but there the be- 
havior of aC/Tc at small and large dopings is generally 
uncorrelated. 

We found reasonably good agreement between 
our theory and the experimental phase diagram of 
Ba(Fci_a;Ruj;)2As2, in which Fe is subsituted by isova- 
Icnt Ru (sec Refsi^ii) and on the data for the doping de- 
pendence of the specific heat jump at T^t^— This agree- 
ment is a good indicator that our theory captures the key 
physics of non-monotonic behavior of AC/Tj,, particu- 
larly the reduction of AC/Tc in the mixed state. Whether 
the data can distinguish between rigid band and disorder 
scenarios is a more subtle issue as the interplay between 
doping-induccd disorder and doping-induced change in 
the band structure is likely to be material-dependent. 
Another subtle issue is the apparent dependence of 
the measured A.C /T,. on both sides of optimal doping. 
Our log-log plots reproduce this dependence reasonably 
well, but our actual AC/Tc are more mild than . One 
possible reason is our neglect of the doping dependence 
of 7 in the normal state specific heat C = 7T. In reality, 
7 also decreases on both sides of optimal doping^ and 



this should sharpen up the T dependence of C^C jT^- 

Finally, there are certainly other elements of the 
physics of Ba(Fei_a;Rua;)2As2 which we neglected in our 
model. In particular, Brouet et al. demonstrated^ that 
Fermi velocities in Ba(Fei_a;Rua;)2As2 are larger than in 
undoped BaFe2As2, this observation likely implies that 
electronic correlations are weaker in Ba(Fci_j;Ru2;)2As2. 
Dhaka et al. argued^ that magnetic dilution due to 
Ru substitution contributes to the destruction of SDW 
order. These effects add on top of Ru-induced impu- 
rity scattering, which we studied here, and call for more 
comprehensive analysis of isolvalent doping of Fe in 122 
materials. 
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